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Abstract 

We study numerically the phase structure of a model of 3D gravity interacting with 
scalar fermions. We measure the 3D counterpart of the "string" susceptibility exponent as 
a function of the inverse Newton coupling a. We show that there are two phases separated 
by a critical point around a c ~ 2. The numerical results support the hypothesis that the 
phase structures of 3D and 2D simplicial gravity are qualitatively similar, the inverse 
Newton coupling in 3D playing the role of the central charge of matter in 2D. 



1 Introduction 



The success of the so-called dynamical triangulations or matrix models as a theory of 
2D quantum gravity and non-critical strings has brought about a hope, that analogous 
discrete approaches might be useful in higher dimensions as well. The most natural way 
to introduce discretized quantum gravity is to consider simplicial complexes instead of 
continuous manifolds. Then, a path integral over metrics can be simply defined as a sum 
over all complexes having some fixed properties. For example, it is natural to restrict 
their topology. In the present paper, we consider only the 3-dimensional case, where, by 
definition, a simplicial complex is a collection of tetrahedra glued along their faces in such 
a way that any two of them can have at most one triangle in common. In addition we 
will restrict the topology to be that of the three-sphere. 

To introduce metric properties, one can assume that complexes represent piece-wise 
linear manifolds constructed by gluing together equilateral tetrahedra [|l| . The volume is 
then proportional to the number of them. In the continuum limit, as this number grows, 
the edge length, a, simultaneously tends to zero: a — > 0. In three dimensions, there is a 
one-to-one correspondence between classes of piece-wise linear and smooth manifolds. In 
other words, any homeomorphism can equally well be approximated by either a piece-wise 
linear or a smooth map. This property makes the model self-consistent. 

Within a piece-wise linear approximation, the curvature is singular: the space is flat 
everywhere except at the edges. Therefore, one should consider integrated quantities, e.g. 
the mean scalar curvature 

( 1 

J d 3 x yfgR = a(2-nN 1 - 6iV 3 arccos -) (1.1) 

where a is the lattice spacing and N , N ly N 2 , N 3 are the numbers of vertices, edges, 
triangles and tetrahedra, respectively. 

For manifolds in 3 dimensions, the Euler character vanishes x — N0—N1+N2—N3 = 0. 
Together with the other constraint N 2 = 2N 3 , it means that only 2 of the iVj's are 
independent, Ni = No + N%, and a natural lattice action depends on 2 dimensionless 
parameters. Now, we are in a position to define the simplicial-gravity partition function 

Z(a,n) = Y.e aNo ~^ = -£Z N3 (a)e-» N > (1.2) 

{C} N 3 

where J2{c} is a sum over some class of complexes (e.g. of a fixed topology, as mentioned 
above), a — 1/G is the inverse Newton constant and /1 is the cosmological constant. 

It can be shown that if N$ is fixed, Nq is restricted from above by iVo < |A^3 + const. 
Therefore, it is natural, keeping a fixed, to let \i approach its critical value, fi c , at which 
the sum over N$ in Eq. ( |1.2|) becomes divergent. In the vicinity of fi c one could hope to 
find critical behavior corresponding to a continuum limit of the model. In general the 
value of fi c will be a function, /x c (a), of the other parameter a. While /i — > fi c (a) in general 
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will correspond to an infinite volume limit (A3 — > 00), it might require a fine-tuning of 
a to find a limit which has interesting continuum physics. From this point of view the 
parameter a might play a role similar to the one played by the (inverse) temperature j3 
in the Ising model. While we can take the infinite volume limit of the Ising model for any 
/3, it is only for one specific value of j3 (corresponding to the critical temperature) that 
we can relate the Ising model to a continuum c = 1/2 conformal field theory. 

For the model to be physically consistent, the free energy per unit volume, log Z Ns (a) , 
has to be finite. It implies that the series over A3 in Eq. (|1.2|) has to have a finite radius 
of convergence: < \i c < +00. If this holds, statistical models defined on such an ensem- 
ble of fluctuating lattices possess reasonable thermodynamical behavior. Thus one can 
hope to describe matter coupled to quantum gravity. This program has been successfully 
carried out in two dimensions with help of the matrix model techniques 0, 0. In the 
3D case, analytic tools are sparse and numerical simulations have so far been the main 
source of information. In the case of pure 3D gravity in references |3|, |J], and in the case 
of matter fields coupled to 3D gravity in references ||. 

Many questions, which are simple in two dimensions, become complicated in three. 
For example, 2-manifolds are classified according to values of the Euler character, and 
therefore topology, homotopy and homology are actually equivalent descriptions. In three 
dimensions, sets of manifolds characterized according to these 3 criteria are essentially 
different. From this point of view it is natural to ask how large a class of complexes {C} 
one can include in Eq. fll.2|) , in order that the model still has acceptable properties, like 
being well defined and having a finite free energy per unit volume. It has been argued 
that it is sufficient to fix only the first homology group H\(C, Z) in order to have a 
finite value of \i c . It can be reformulated as the following statement: "The number of 
simplicial manifolds constructed of the given number of tetrahedra, N 3 , and having a fixed 
homology type (viz., the first Betti number and the torsion coefficients) grows at most 
exponentially function of A/3" . 

Let us further note that the existence of such an exponential bound will remain valid 
if one attribute positive weights which grow at most exponentially with A/3, to every 
complex. The other way around, if one manages to represent a sum over complexes as 
a weighted sum over another class of objects, the weights again exponentially bounded, 
then it is sufficient to prove that the number of objects in the new class is exponentially 
bounded. In the class of regular triangulations there is firm numerical evidence for an 



exponential bound ||10|| . Further, analytical arguments now exist, which strongly favor an 
exponential bound |TT] . 

The simplest example of a weighted sum of complexes is given by the canonical parti- 
tion function Zn 3 (ci) introduced in Eq. ( |1.2| ). Another example of physical interest is the 
partition function in the presence of free matter fields. To introduce it let us consider the 
set of polyhedral complexes dual to simplicial ones. Their 1-skeletons are some 4-valent 
graphs, whose vertices correspond to tetrahedra, and links are dual to triangles. 
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Such a graph can be defined by the adjacency matrix 

^.(4) f 1 if vertices i and ?' are connected by a link .„ . 

G ii = \ n 4.1. • L3 

J [ U otherwise 

Of course the set of 4-valent graphs is not identical to the totality of simplicial com- 
plexes. Similarly, in two dimensions, the set of ordinary <f) 3 Feynmann diagrams is different 
from that of triangulations - one has to introduce the notion of fat graphs to establish 
the equivalence. 

If one attaches a n-component vector xf (fi — 1, . . . , n) to the z-th vertex, the corre- 
sponding gaussian integral can be performed explicitly for every given graph G^: 

Z^UG^) = J II II^e x P[-i E Gf«-^) 2 ] = (det'L«) *, (1.4) 

where the discrete Laplacian is given by 

L$ = A6 tj - Gf . (1.5) 

There is the nice combinatorial representation of the determinant given by Kirchoff's 
theorem [|15|1 : 

det'L = \T(G)\, (1.6) 

where T(G) is the set of all possible spanning trees embedded into a graph G or, equiva- 
lently, the number of connected trees which can be obtained from the graph by cutting its 
links. To remove the zero mode in Eq. ( |1.4|) , we fixed the field at the iVyth vertex. Kir- 
choff's theorem states that det 'L does not depend on the choice of the vertex. Therefore, 
|T(G)| is the number of rooted trees. 

The number of spanning trees of a fc-valent graph G^' with n vertices can be estimated 
from above as |15j 

\T(G^)\ < -f^r 1 . (1.7) 
n \n — 1 J 

The case we are particularly interested in is the 2-component Grassmann field (or 
equivalently, n = — 2 in eq. (|1.4| )), where we obtain the gravity+matter partition function 
in the form 

Z«(a,/i) = ^det'Le^"^ 3 = ^zgHe^. (1.8) 

{C} N 3 

The weights, det 'L, are positive integers, hence, 

Z Na (a) < Zjg(a). (1.9) 

In two dimensions, this type of matter has the central charge c = —2 and the corre- 
sponding matrix model has been solved explicitly by purely combinatorial means || . We 
can repeat the same trick in 3 dimensions: 

ZW(a,/i) = ^Tdet'Le^-^ 3 = E E e aNo ^ N:i = E E e aNo ~^ Ns . (1.10) 

{°} {°} {T(C)} {?} {C(T)} 



3 



Here we have a sum over complexes, D{c}> weighted with the number of spanning trees 
in the dual 1-skeleton of each complex, J2{f(c)}- ^ e can obtain exactly the same class of 
configurations by taking the sum over all possible trees, X^-p weighted with the number 

of complexes, Ylsc(T)y wmcn can be reconstructed from a given tree T by "restoring" links 
which we imagine have been cut in order to produce the tree from the original graph. 

In three dimensions, the trees appearing this way correspond to spherical simplicial 
balls obtained starting with a single tetrahedron and subsequently gluing other tetrahedra 
to the faces of the boundary If we denote by n , ri\ and n 2 the numbers of vertices, edges 
and triangles on the boundary of a ball, we find that 

n = N 3 + 3, n 2 = 2(N 3 + l). (1.11) 



It is convenient to weight every complex with the number of spanning trees T. This is 
an exponential reweighting of the kind discussed above. It can be done by introducing an- 
other system of free matter fields attached to vertices of triangulations. Then, analogously 
to Eq. ( 1.10|) , we define the modified partition function 



Z( 2 )(a,/i) = ]T det 'L det 'L e aN ^ Na . (1.12) 

{C} 

This partition function was used in || to prove that the number of simplicial 3-manifolds 
having a fixed homology type grows exponentially with the number of tetrahedra, as men- 
tioned above. A priori it is not clear if Z^^a, /1) defines a theory in the same universality 
class as the one defined by Z^(a,fx). However, Z^(a,fx) is the partition function of a 
set of Gaussian fields coupled to geometry in a well-defined way, and it is from this point 
of view equally interesting as the system defined by ZW(a,/i). In the next section we 
consider a reduced version of ( |1.12j ). 



2 Reduced model and connection with the 2-dimensional 
loop gas model 

In 3 dimensions there are several natural choices of classes of complexes one can use in the 
definition of the partition function, as already discussed, and the question of universality 
classes is open. A simple choice is the set of 3-spheres with imbedded so-called 2-collapsible 
spines. A spine is the 2-skeleton, K 2 , of a complex C in a cell-decomposition with one 
3-cell. In the simplicial context this 3-cell is a ball -B3 with a triangulated boundary. Thus 
we are lead to study the following partition function: 

zS(«,aO = E E E , (2-i) 

C~S 3 K 2 CC TcK 2 

where J2c~s 3 * s the sum over all simplicial 3-spheres, where J2k 2 cc * s the sum over 2- 
collapsible spines, K 2 , embedded in C, and where J2tck 2 * s the sum over all spanning trees 



4 



in K 2 . The terms in fl2.1|) is a subset of the terms appearing in (|1.12|) when rewriting the 
determinant using a trick similar to the one used in ( |1.10| ) (see M for details). The spine 
K 2 in ( p.l| ) is obtained from the 3-sphere C a pair-wise identification (gluing) of all the 
boundary triangles. If K 2 is 2-collapsible, the pair-wise identification can be made "local" 
such that one can use recursively the gluing operation consisting of the identification of 
2 triangles sharing a common link on the boundary, or, in other words, a folding along 
the link. In the dual language, one starts with an arbitrary spherical 3-valent fat graph 
and applies the move which can be represented as the flip of a link with the subsequent 
elimination of it: 



(2.2) 

Both of these moves are well known and have been used in Monte-Carlo simulations of 
triangulated surfaces 0. Having glued all triangles, one finishes with a collection of 
self-avoiding closed loops, the number of which equals N — 1. 

It is convenient, instead of erasing a link after a flip, to decorate it with a dashed 
propagator and keep track of it in the process of subsequent foldings. To represent such 
configurations, we need to introduce the infinite set of vertices 



and so on 



(2.3) 



which can be generated by successive flips. For example, 



gives 



gives 



or 



or 



give 




give 



and so on. 

In the end, we obtain planar diagrams with two types of propagators: solid ones 
produce closed loops while dashed form arbitrary clusters of dotted lines. The total 
number of vertices is equal to 2(N 3 + 1). We weight every loop with the numerical factor 
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Thus, the model can be represented as a non-local gas of closed loops with the partition 
function 

ZS(«,A0 = £^^-^, (2.4) 

{V} 

where Yl{v} ls the sum over all the diagrams. The weight w-p is equal to the number of 
inequivalent starting configurations giving the same planar diagram T> (can be 0). The 
number of closed loops is simply equal to No — 1 and a can be identified with the inverse 
Newton constant. 



We see that the weights wp in Eq. Q2.4|) in principle are recursively calculable. How- 
ever, from a practical point of view it is difficult to take into account the entropy of the 
configurations, and the model Q2.4J) seems to be analyticlly unsolvable. Therefore, let us 
consider the localized version of the model, namely, the dense phase of the self-avoiding 
loop gas matrix model (in the following simply denoted the loop gas (LG) model): 

AT _ N " „„ N 



f /V /V /V 

7j ~lg = /^Vn^%exp[--trF 2 -^ytrX, 2 + ^-ytrFX, 2 

J v=l A A v=l 1 u=l 

= yd w2 yexp[-ytry 2 -^tr 2 log(l®l-Aty®l-l®Aty)] . (2.5) 



Here Y and X v are hermitian N x N matrices. We have attached the lower index to the 
gaussian X v variable to weight every closed loop with the factor n = e a . In the N — > oo 
limit, this model generates planar diagram having only the simplest 3-valent vertices from 
the whole series (|2.3| ). For such diagrams, all the weights in Eq. (|2.4|) are trivial: w-p = 1. 

This truncated local model might have some qualitative features in common with 
3-dimensional simplicial gravity. It seems to be a reasonable assumption, because the 
corresponding universality class is very large, (see Ref. 0), and includes all local per- 
turbations of the model ( |2.5| ). On the other hand, multiple overlapping sequences of the 
flips can be regarded as a kind of renormalization group procedure, which could draw any 
given system to a stable fixed point in the space of all planar-graph models. It looks a 
bit more natural in dual terms. One takes an arbitrary spherical ball with a huge tri- 
angulated boundary and uses subsequently and randomly the folding operation (dual to 
the flip). Soon the boundary triangulation will be randomized. If this randomization has 
appropriate robust statistical properties, we would expect to find many models falling in 
the same universality class. 

Then any 2 models from the same universality class should be identical in their con- 
tinuum limits. 

The critical behavior of the loop gas matrix model is well known. For n < 2, it 
describes 2D gravity interacting with c < 1 conformal matter 

c = 1 — 6- — — ; n = — 2cos7rg , (2.6) 

while for n > 2 the corresponding matter is non-critical and the model trivializes except 
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in special situations fL2f . In this phase, the number of closed loops is proportional to the 
volume and the mean length of each remains finite. 

In the vicinity of a critical point u c , the partition function behaves as 

Z LG « Gu c -/i) 2 -^. (2.7) 

This formula defines the famous string susceptibility exponent 7 str . 
Two main quantities of interest are 

(N 3 ) LG = |-logZ LG « _( 2 - 7str )-^- (2.8) 

and 

(N ) LG = ^-\ogZ LG « (2- 7str ) f c(w) - log |^- Ai| 

« --(iV 3 )LG + ^log(iV3) LG (2.9) 

where the coefficients can be calculated explicitly using the Gaudin and Kostov's exact 
solution [0] 

d 11 d% tI 1 

— logHc = -- — — ; — — = - — — r. (2.10) 

on 2 n + 2 an sin 7r(l — g>oJ 

If 2 - n < 1, then 

(N )lg ~ 7(^ 3 )lg + i log(iV 3 ) LG . (2.11) 

4 7rv2 — n 

If n > 2, then 7 str is independent of n and (N Q ) lg is proportional to (N 3 ) lg, except in 



special situations [12 



3 Numerical experiments 

Unfortunately, the nonlocal model ( p.4| ) looks analytically unsolvable. However, it can be 
successfully simulated numerically by the recursive sampling-Monte-Carlo technique. It 
is natural to use a recursive sampling to perform the summation over the tree-dual balls 
in a close analogy with Ref. [T(J. Then, a Metropolis algorithm can be used to produce 
a random walk in a space of all simplicial spherical balls. Closed manifolds appear when 
the boundary of a ball is shrunk, again by a trial and error Monte Carlo procedure. In 
the following we describe the algorithm in some detail. 



3.1 The algorithm for generating three-spheres 

As mentioned, a recursive algorithm has been used with great success for a c = —2 model 
coupled to two-dimensional quantum gravity. Since our model resembles the c = — 2 
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Figure 1: Graphical presentation of Schwinger- Dyson equation for 3-valent tree-diagrams. 

model quite a lot, it is natural to try to use the same technique in three dimensions. 
One of the main reasons for the success of the algorithm in the two-dimensional case 
was that the branching probability of tree-graphs and the counting of rainbow diagrams 
could be calculated analytically. In the present three-dimensional model the branching 
probability can be calculated analytically as shown below. However, while there is a 
three-dimensional analogy of the rainbow diagram encountered in two dimensions, they 
are not two-dimensional extended objects and we do not know how to count them. 

Let us calculate the branching probability of 3-valent tree-diagrams. The Schwinger- 
Dyson equation for 3-valent graphs corresponding to Fig. 1 is 

T = xT 3 + l, (3.1) 

where T(x) is the generating function for the number of tree diagrams, i.e., 

T(x) = £t n x n , 

n 

where t n is the number of 3-valent tree-diagrams with n vertices. 

Since the Schwinger-Dyson equation is a cubic equation with respect to T, it has three 
solutions. We have to select the unique solution which satisfies the physically imposed 
boundary condition t — 1. This results in the following solution: 



T = *mhh(±]og(z + y/T+*j), * = iffi, (3-2) 
which has the following power series expansion: 

T(x) = Vt/ = V ,/ 1 - x n . (3.3) 

K> to ^n!(2n + l)! 1 ] 

We can now define the probability of the 3-valent tree-diagrams with k vertices branch- 
ing into tree-diagrams with /, m and n vertices 

w(l,m,n) = 1 m - (l + m + n — k — 1), (3.4) 
tk 

which then naturally satisfies 

y^w(l, m, n) = 1, (3.5) 
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where the summation is over l,m,n > with the constraint l + m + n — k — 1. By 
rearranging the summation we obtain the following relation 

w(l,m,n) = Wi(l) W2(m, n), (3.6) 

where 

/ a w(l,m : n) 

= w(l,m,n), w 2 {m,n) = ryr — , (3.7) 

m',n'>0 Wl ") 
m' + n' = fe - I - 1 

where m + n — k — I — 1. Since we have the relations 

fc— i 

5>i(Z) = 1, E ^ 2 (m,n) = 1, (3.8) 

1=0 m,n>0 

m + n = k — I — 1 

we can divide the branching process into two parts using the relation 

E ™(Z,m,n) = 3( V-Z~ 2 (O^^^ 1 )' (3-9) 

m, n > ' 
m + n = fc — I — 1 

and write: 

3(/c — /) — 2 titk-i-i k — l t m t n . . 

^i(0 = 7—^ — , w 2 (m,n) = — . (3.10) 

Using these analytic formulas we can set up a recursive sampling algorithm for construct- 
ing 3-valent tree-diagrams in three dimensions. 

In order to construct a three-dimensional simplicial manifold with the topology of S* 3 
from the 3-valent diagrams, we need to identify the "end branches" of the tree-diagrams 
in the dual diagram sense, or equivalently, successive neighboring triangles of the corre- 
sponding original manifold. In contrast to the two-dimensional c = —2 model, where an 
analytic formula for the number of rainbow diagrams is given, it is presently unknown how 
to obtain a closed expression in the three-dimensional case since the trees are now gen- 
uine three-dimensional objects. The initial three-dimensional tree diagram corresponds 
to the tree shape volume object which has a ball topology with a triangulated S 2 sphere 
as a boundary. It should be noted that in this model all the links of the corresponding 
three-dimensional simplicial manifold are on the boundary in the initial tree configura- 
tion. We first randomly pick up a link on the boundary and identify the neighboring 
triangles which share the link as a common one. This link is now hidden inside from the 
boundary. We then pick up another link which is on the boundary surface and repeat 
the gluing procedure of the triangles. This procedure should be continued successively 
until all the triangles on the boundary surface are identified (and the boundary thereby 
has vanished). Since we don't know the analytic expression of the weight factor of the 
three-dimensional "rainbow diagrams", we need to introduce the ungluing process and 
use a Monte Carlo simulation to obtain a typically triangulated three-dimensional S 3 
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sphere in our model. Since the gluing process is carried out by a Monte Carlo algorithm, 
we operate both with the gluing process and its reciprocal in such a way that the usual 
detailed balance is fulfilled. In this procedure the order of gluing and ungluing processes, 
which are reflected in the internal configuration of the simplicial manifold, are important 
and thus should be memorized in the simulation. In this Monte Carlo procedure there is 
unavoidably introduced an iteration uncertainty in the counting procedure. 



3.2 Observables 

After constructing typical S 3 manifolds one can use these to calculate expectation values 
of observables. A basic observable which probes the phase diagram of the model is the 
"string" susceptibility exponent 7 defined analogously to Eq. ( |2.7| ). The most convenient 
way to extract 7 is by means of the so-called baby universe technique ]17| , which contrary 



to the original definition (|2.7|) allows us to extract 7 from an ensemble of manifolds with 
fixed volume A3. One counts the configurations of three-manifolds which are separated 
into two parts via a common triangle, and then measures the volumes M 3 and A3 — M3 of 
the separated parts. The string susceptibility is then derived from the following relation: 



W(M , M 3 Z{M 3 ) (A3 - M 3 )Z{N 3 - M 3 ) 
W(M 3 ) = N 3 



■M 3 / 


M 3 \~ 


7-2 


.N 3 y 


1 -) 

a 3 K 





(3.11) 



where the partition function Z(N 3 ) (in accordance with ( |2.7| )) is assumed to be parametrized 
as Z(N 3 ) ~ A^3~ 3 e~^ cNi . In particular we obtain the following approximate scaling func- 
tion: 

= W»£^~( 7 -2)ri--!-V (3-12) 



W(M 3 ) v ' ' \x 1 

where the scaling parameter is x = M 3 /N 3 . 

A number of measured values of 7 as a function of a, the inverse Newton constant, 
is shown in Fig. pj. Statistical error-bars are shown in the Figure. The two sub-figures in 
Fig. |2| represent independent measurements for the different value of N 3 and they agree 
perfectly within error-bars. Systematic errors {e.g. due to finite size effects) seem to be of 
the same order except in a vicinity of the transition point. The data are collections of three 
simulations which took place approximately 4-6 months of CPU time for each calculation 
on a Pentium 200 MHz personal computer and HP Convex and Origin 2000 parallel 
computers. We have collected data for A 3 = 50, 100, 200 at log(l/2) < a < log(256), 
and A 3 = 300 at a = and have shown the extracted 7 in the upper sub-figure of Fig. ||. 
Other data for A3 = 300, 500 lead the 7- values shown in the lower sub-figure. Although 
the number of tetrahedra are small, they seem to give reliable measurements of 7 in the 
interval —0.7 < a < 5.5, except for A3 = 50. For the value of A3 there is a systematic 
deviation from the other values, probably caused by the small number of A3. A typical 
distribution of logiy(M3) for a = and A 3 = 300 is shown in Fig. [3] where the slope of 
the curve determines 7. 
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Figure 2: "String" susceptibility 7 versus a, the inverse Newton constant, for N3 = 50(o) 
100(-), 200(D), 300(x) above and N 3 = 300(o), 500(-) below. 
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Figure 3: log W(M 3 ) versus log (M 3 (N S - M 3 )) at a = for iV 3 = 300. 
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Figure 5: log-^j^^ versus M 3 at a 



for iV 3 =300. 



The scaling function F(x) with the predicted curve (|3.12| ) is shown in Fig. |] for a = 
with N 3 = 50, 100, 200, 300. We can see that all the data scale nicely within the error- 
bars except for M 3 <^ M™ m = 5. As the formula (|3 . 1 2|) is valid only asymptotically, we 
can use fitting only for values of M3 bigger than some cutoff M™ m . The independence of 
results from this cutoff can serve as an indicator of the reliability. A corresponding typical 
measurements are shown in Fig. For small values of M3 < M™ m = 5, the measured 
"virtual" values of 7 oscillate and for larger values they gradually converge to a stable 
limit. 



In the measurements of 7 we have thus collected the data satisfying M3 > M™ m +l = 6. 
We have measured 7 with good statistics for a = where we collected 10 5 configurations 
while for the other values we collected about 2 x 10 4 configurations. The measured value 
7 is so close to —1 that it is tempting to conjecture that this is the exact value for a = 0. 
One could hope that it was possible to prove this result analytically, since it corresponds 
to the strong coupling limit G = 1/a = 00. 



4 Discussion 

The results of the numerical experiments give some support to the idea that there should 
be a qualitative similarity of the phase structures of 2- and 3-dimensional simplicial gravi- 
ties with Newton's coupling constant playing a role analogous to the matter central charge 
in 2 dimensions. A similar idea can be found in a slightly different context in |TJJ. The 



experimental values of 7 shown in Fig. |2| support the picture of 2 phases with the critical 
coupling « c ^ 2. For a < a c , the simplest square root function 



7 



VI -0.64a (4.1) 
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fits surprisingly well the experimental points. We show the above expression by the dotted 
curve in Fig. |2|. For a > a c , 7 is a constant within error bars about 0.3. Since the systems 
are quite small, and since there still is some drift in the extracted values of 7 towards 
larger values when N 3 is increased (see in particular lower sub- figure in Fig. |2|), we cannot 
conclude if 7 = 1/2 in the limit of infinite N 3 or if 7 = 1/3, which is the next lowest value 



that 7 can have in the branched polymer phase |I3 . 

From a theoretical point of view, it has been recognized by H. Kawai and one of 
authors (Y.W.) that in the case of D > 2 and around I /a = G = 0, the branched 
configurations give much more dominant contribution to the path-integral over metrics 



than the configurations near the flat space-time ||18||. Let us consider the scaling property 
of 



X = expj^- jy/gR). (4.2) 

Near flat space-time dimensional analysis is reliable and leads to Jy/gR ~ V 1 ^^ . Then 
one finds 

( V l -T>\ 

X ~ expl const. — — — J , (4.3) 

where V is the volume of space-time. On the other hand, for branched polymer like 
configurations it is easy to show that i.e. Jy/gR ~ V. Thus one finds in this case that 

X ~ exp ( const.— ^ (4.4) 



instead of (|4.3|) . Comparing eqs. (|4.3|) and (|4.4|) around G ~ 0, one finds that (|4.4j) gives 



a more dominant contribution to the path-integral than the contribution ( f4.3|) from flat 
space. This is one reason why one cannot perform the quantization of gravity perturba- 
tively around the flat space-time (except for the very special cases when l/a c — G c — 0). 
This theoretical observation is consistent with the result of this paper. 

If we accept the result that there are two phases in the present model of 3-dimensional 
simplicial gravity, then the next natural question to ask is: "How large is the corresponding 
universality class?" Does the complete model ( |L12j ) with a sum over all 3-spheres belong 



to the same class? It seems plausible, but it is difficult to judge how restricted the subset 
of simply connected 3-manifolds taken into account in the simulated statistical ensemble 
really is. It is not easy to construct a simplicial 3-sphere which would not be included 
in it. At first thought, Zeeman's Dunce Hat[19j could serve as a mean to construct a 



counter example. Triangles belonging to a boundary of an initial ball cannot form such a 
pseudosurface after local foldings. However, those lying inside the ball can. Even in the 
simplest case of 2 tetrahedra, a Zeeman's Dunce Hat can easily be obtained. Further, the 
construction given in Ref. is not a "counter example" because it is not a triangulation. 
On the other hand, there is an old result of Haken |2l| who showed that any 3-sphere 



possesses an obviously trivial presentation of tt\ (in the sense that it can be trivialized by 
a sequence of simple substitutions of the length less or equal than 4). Unfortunately, this 
result does not imply that any simplicial 3-sphere allows for at least one 2-collapsible spine 
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(if it did, the famous Poincare hypothesis would be simultaneously proven). It seems an 
open question whether any simplicial 3-sphere is counted at least once in the partition 
function ( |2.4| ). If it is the case, then only the matter sector of the original model ( |1.12| ) 
is distorted by our simplifications leading to the reduced model. Even if not all 3-spheres 
are included in the model, it seems reasonable to conjecture that the reduced model and 
the original model belong to the same universality class as far as their continuum limits 
are concerned. 
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